Stark-tuned Forster resonance and dipole blockade for two to five cold Rydberg 
atoms: Monte Carlo simulations for various spatial configurations 
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Results of numerical Monte Carlo simulations for the Stark-tuned Forster resonance and dipole 
blockade between two to five cold rubidium Rydberg atoms in various spatial configurations are 
presented. The effects of the atoms' spatial uncertainties on the resonance amplitude and spectra 
are investigated. The feasibility of observing coherent Rabi-like population oscillations at a Forster 
resonance between two cold Rydberg atoms is analyzed. Spectra and the fidelity of the Rydberg 
dipole blockade are calculated for various experimental conditions, including nonzero detuning from 
the Forster resonance and finite laser line width. The results are discussed in the context of quantum- 
information processing with Rydberg atoms. 

PACS numbers: 32.80.Ee, 03.67.Lx, 34.10.+X, 32.70.Jz 



I. INTRODUCTION 

Strong dipole-dipole interaction (DDI) between highly 
excited Rydberg atoms with a principal quantum num- 
ber n 3> 1 is a key element for quantum logic gates with 
qubits based on single alkali-metal atoms [1,2]. Two ba- 
sic proposals consider a short-term DDI of two close Ry- 
dberg atoms to implement conditional quantum phase 
gates [3] or laser excitation of only one Rydberg atom 
in a mesoscopic ensemble (dipole blockade) to obtain en- 
tangled qubits [4]. 

DDI arises as soon as two ground-state atoms are laser 
excited to the Rydberg states of the opposite parity [2,5]. 
Alternatively, two Rydberg atoms in an identical state 
may interact via a Forster resonance if the atoms are ex- 
cited to a level that lies midway between two other levels 
of the opposite parity [6] . Due to DDI one of the atoms 
undergoes a transition to a higher Rydberg state, while 
the other atom goes to a lower Rydberg state. This pro- 
cess is resonantly enhanced at zero energy defect. For 
some Rydberg levels Forster resonances can be precisely 
tuned with an electric field via the Stark effect [2,6]. Ac- 
cidental quasi-Forster resonances with a small energy de- 
fect can also be found for certain Rydberg states [7,8]. 

Accidental quasi-Forster resonances are most suitable 
for the dipole blockade, which is insensitive to decoher- 
ence and fluctuations of the DDI strength [1,8]. Dipole 
blockade directly generates entanglement between qubits. 
Two recent experiments [9,10] have demonstrated entan- 
glement and controlled-NOT quantum gates for the two 
optically trapped and individually addressed Rb atoms. 
Both experiments employed the dipole blockade at the 
laser excitation to the Rydberg states interacting via 
quasi-Forster resonances. The energy defects of a few 
MHz were small compared to the DDI strength, and 
blockade efficiency was close to 1. Partial dipole blockade 
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has also been demonstrated in several earlier experiments 
with quasi-Forster or even van der Waals interactions be- 
tween Rydberg atoms [11-14]. 

An alternative approach to quantum-information pro- 
cessing with Rydberg atoms relies on the conditional 
quantum phase gates (QPG) [7,15,16]. This approach 
has the advantage that it does not require individual ad- 
dressing of the atoms and strong DDI [1], but it requires 
coherence to control the phase of the collective wave func- 
tion of the two atoms. Coherent interaction means that 
Rabi-like population oscillations of high contrast between 
the initial and final collective states must be obtained. 
This can be achieved with the resonant DDI, in partic- 
ular, with Stark-tuned Forster resonances. Compared 
to the accidental quasi-Forster resonances, Stark-tuned 
resonance is more flexible in controlling the interaction 
strength. A weak electric field can be used for fast switch- 
ing of the interaction [17] and laser excitation [18], and 
for adjusting the phase of the collective wave function by 
inducing the Rabi-like population oscillations and stop- 
ping them at a proper time. Efficient dipole blockade can 
be implemented with Stark-tuned Forster resonances, as 
well. 

QPG or complete dipole blockade at a Stark-tuned 
Forster resonance between two or a few Rydberg atoms 
has not been reported yet. For many atoms, partial 
dipole blockade was observed [19,20], and experiments 
on the population dynamics at Forster resonance have 
been done [21,22]. 

In our recent experiment [23] we observed for the first 
time a Stark-tuned Forster resonance between two rubid- 
ium Rydberg atoms randomly positioned in a small laser 
excitation volume. High-resolution spectra of the Forster 
resonance have been obtained for N = 2 — 5 of the de- 
tected Rydberg atoms that allowed us to investigate in 
detail the line shape dependence on N and effective co- 
herence time. The related theoretical analysis relied on 
numerical Monte Carlo simulations for two to five Ryd- 
berg atoms. Although the theoretical model was simpli- 
fied by ignoring the hyperfine structure of the Rydberg 
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levels and the broadenings due to parasitic electromag- 
netic fields, we have shown that these broadenings can 
be taken into account by reducing the effective interac- 
tion time. This reduction results in a Fourier broadening 
of the resonance to a value that matches the experimen- 
tal width. A good agreement between experiment and 
theory has confirmed the validity of the model used. 

In this article we present a more detailed description 
of our Monte Carlo model and apply it to simulate a 
possible observation of the coherent Rabi-like population 
oscillations at a Stark-tuned Forster resonance. We also 
simulate numerically the excitation spectra for the dipole 
blockade effect at a Forster resonance. The simulations 
were performed for various spatial configurations of the 
interacting Rydberg atoms relevant to the typical exper- 
imental conditions. The effects of the atoms' spatial un- 
certainties on the Forster resonance and dipole blockade 
are analyzed. 



II. INTERACTION OF TWO RYDBERG ATOMS 

A. Two frozen atoms 

The operator of the dipole-dipole interaction between 
two atoms a and b is 
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FIG. 1: (Color online) (a) Scheme of the relevant energy lev- 
els and transitions in Rb Rydberg atoms. Laser radiation 
initially excites the 37P state from the 55 ground state. Res- 
onant dipole-dipole interaction (Forster resonance) induces 
transitions to the neighboring 375 and 385 states, (b)-(d) 
Typical spatial atom distributions in experiments with two 
interacting Rydberg atoms: (b) two atoms randomly placed 
in a single interaction volume, (c) two atoms in the individual 
cigar-shaped dipole traps, and (d) two atoms in the individual 
dipole traps evenly localized in all directions. 



state. In this case a single Forster resonance is observed 
at 1.79 V/cm [23]. Single resonance is most appropriate 
for the experimental and theoretical analysis and for im- 
plementing the quantum gates, because the overlapping 
of several resonances leads to decoherence. 

Our Forster resonance induces transitions between Ry- 
dberg states with AMj = 0. This corresponds to the z- 
oriented dipoles, and Eq. (fTJ) reduces to a simpler form: 
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Here d Q and d b are the dipole-moment operators for the 
atoms a and 6, R a b is a vector connecting the two atoms, 
and £o is the dielectric constant. The interaction en- 
ergy depends on the mutual orientation of the two atom 
dipoles (defined by their quantum states) and on their 
orientation with respect to R a &. 

Degenerate atomic states usually produce several 
Stark-tuned Forster resonances corresponding to dif- 
ferent momentum projections. Nevertheless, various 
Stark-tuned Forster resonances in Rydberg atoms should 
have similar spectroscopic properties. The difference is 
mainly in the dipole moments that define the interaction 
strength and in the initial resonance detuning in a zero 
electric field. Therefore, the calculations can be done 
for an arbitrary Forster resonance, whose theoretical de- 
scription can be scaled later to other dipole moments and 
initial detunings. 

In this article we consider an example of the Stark- 
tuned Forster resonance Rb(37P 3 / 2 ) + Rb(37P 3 / 2 ) — > 
Rb(375 1/2 ) + Rb(385 1/2 ) for two or more Rb Ryd- 
berg atoms [see Fig. 1(a) for the level and transition 
scheme]. The initial energy detuning [2E(37P 3 / 2 ) — 
E(375i /2 ) - E(385i /2 )]//i in a zero electric field is 103 
MHz. The proper choice of the exciting laser polariza- 
tion (along the dc electric field) provides excitation of 
only 37P 3 / 2 (|Mj|=l/2) atoms from an intermediate 65 
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where d a ^ are the z components of the dipole-moment 
operators of the two atoms and Z ab is the z component 
of R a b (z axis is chosen along the dc electric field). 

For the short interaction time to < 1 fis, which is of 
interest for fast quantum gates, we can ignore the sponta- 
neous decay and blackbody-radiation-induced transitions 
from our Rydberg states (the effective lifetimes are tens 
of microseconds for n > 30 [24] ) . This also simplifies the 
calculations by ignoring the motion of the atoms (dis- 
placement < 1 /zm for atom temperatures T < 200 /iK) 
and the hyperfine structure (< 500 kHz for n > 30). 
Then, in order to calculate the time evolution of the 
two Rydberg atoms, we can simply solve the Schrddinger 
equation for a quasimolecule formed by two atoms. For 
the motionless atoms this problem is solved analytically. 
Let us denote the lower energy state 375 as state |1), the 
middle state 37P 3 / 2 as state |2), and the upper state 385 
as state |3), as shown in Fig. 1(a). The collective wave 
function for two Rydberg atoms is 



# = a 22 {t)\22) + a 13 (i)|13) + a 3 i(t)|31), 



(3) 



where 1 22) is the initial collective state populated by 
a short laser pulse at t = 0, 1 13) and |31) are the 
two equally populated final states having a small en- 
ergy detuning A from state 1 22) , and a.y (i) are the 
time-dependent probability amplitudes to be found. The 
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ground state 5S, denoted as state |0) in Fig. 1(a), is ig- 
nored here as it does not contribute to the Forster reso- 
nance. We also ignore the other collective Rydberg states 
|21), |12), 1 23), 1 32} . These states have large energy de- 
tunings from state 1 22) and are not populated by DDI at 
the Forster resonance. A short laser pulse implies that 
the Fourier width of the laser radiation is large and it 
precludes the dipole blockade effect at laser excitation. 

The signal measured in our experiments is a fraction of 
Rydberg atoms in final state |1) or a population of final 
state |1) per atom, which is calculated for the interaction 
time to as 



P2 = ^(k3(io)| 2 + |a 3 i(io)| 2 ) 



(4) 



The solution of the Schrodinger equation with the initial 
conditions 022(0) = 1 and 013(0) = 031(0) = then yields 



atoms into close dipole traps evenly localized in all direc- 
tions [Fig. 1(d)], like in two- or three-dimensional optical 
or magnetic lattices at low atom temperatures (see, e.g., 
Ref. [25]). These configurations appear to be most typi- 
cal for the future experiments on quantum-information 
processing. We are interested in analyzing how their 
intrinsic atom position uncertainties would modify the 
spectrum of the ideal Forster resonance given by Eq. ([SJ . 

No analytical solution can be obtained for the arbi- 
trary spatial distributions. However, we have found that 
for a disordered ensemble of Fig. 1(b) it is possible to ob- 
tain a formula for the spatially averaged resonance am- 
plitude (|02(A = 0)). In order to do this we need to 
choose a proper spatial distribution of the atoms to de- 
fine probability P{r) to find two atoms at distance r in 
the excitation volume. The nearest-neighbor probability 
distribution by Chandrasekhar [26] is appropriate for two 
atoms: 
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where Q a b = V a b/h is the DDI energy in the circular 
frequency units, and A = (2E2 — Ex — E3) /H is the de- 
tuning from exact Forster resonance, which is varied by 
the electric field. The interaction between two motion- 
less Rydberg atoms thus leads to population oscillations 
in the quasimolecule similarly to the Rabi oscillations in 
a two-level atom. Since this process is coherent, it can 
be used to implement two-qubit QPG. At A = the 
first oscillation minimum in Eq. ([5]) corresponds to the it 
phase shift of the collective wave function, which is neces- 
sary for QPG. Experimental observation and study of the 
Rabi-like oscillations is thus a prerequisite for quantum- 
information processing with neutral atoms. 



B. Two spatially fluctuating atoms 

In the realistic experiments the atom positions are not 
fixed even for the cold atoms in tightly focused dipole 
traps, due to finite atom temperatures [16]. The inter- 
action strength is therefore a fluctuating value that may 
result in decoherence and washing out the Rabi-like os- 
cillations in Eq. ([5]). To calculate the signals measured 
in the experiments, Eq. ([5]) must be averaged over all 
possible spatial positions of the two interacting atoms. 

In the experiments with cold trapped atoms there are 
three spatial distributions that deserve special consider- 
ations. The first distribution is a completely disordered 
ensemble of the atoms confined in a small single exci- 
tation volume [Fig. 1(b)], which is formed by two in- 
tersecting laser beams, as in our experiment [23]. The 
second distribution is trapping of the individual atoms 
into close cigar-shaped dipole traps obtained at the fo- 
cal points of the two laser beams [Fig. 1(c)], as in Refs. 
[9,10] . The third distribution is trapping of the individual 
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where r$ ss [3/(47rno)] is the average distance between 
nearest- neighbor atoms at volume density uq. It is valid 
for the atoms homogeneously distributed in the excita- 
tion volume. 

Another simplification we use is a conventional removal 
of the orientation-dependent part in Eq. ([5]); that is, we 
consider a " scalar" dipole-dipole interaction. We remove 
the 3Z 2 b /i?^ b term in Eq. ((2]), so that the interaction 
strength is now given by 1 / ' B? a b- Its validity for disordered 
atom ensembles is grounded by the fact that mutual atom 
orientations are effectively averaged in such ensembles. It 
will be tested later by comparing with numerical results. 
Then the resonance amplitude in Eq. (JS|) is averaged as 
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t a P(r)dr, (7) 



where d2i and d23 are the dipole moments of transitions 
1 2) — > |1) and |2) — > |3). For convenience, we replace the 
expression in the brackets with Oor^/r 3 , where 



\Z2d 2 lfl'23 i 
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is the interaction pulse area at the average distance ro- 
Further simplification of Eq. (|7|) is made by the replace- 
ment x = r 3 /rg that yields 



WA = 0)} 
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An approximate analytical solution may be found if we 
note that the main contribution to the integral in Eq. Q 
comes from the very short distances, where rapidly oscil- 
lating squared sinus function can be replaced with its 



4 



average value 1/2. Then the integration should start at 
x = and stop at some point x c « af9 , where a ~ 1, be- 
cause the integrand exponentially drops at x > x c . This 
approach finally gives a simple result: 



(p 2 (A = 0)) w - (l - e" 



(10) 



This is a general formula that describes both weak dipolc- 
dipole interaction (amplitude grows as a#o/4) and strong 
interaction (amplitude saturates at 1/4) . It may be used 
to calculate resonance amplitude for the arbitrary dipole 
moments, atom densities and interaction times, which are 
all included in 9q. The only fitting parameter is a, which 
can be found by comparing Eq. (fTU)) with more precise 
numerical calculations. 



C. Monte Carlo averaging for two disordered atoms 

Validity of Eq. (fTU|) has been tested by numerical 
Monte Carlo simulations. In this approach, the two Ry- 
dberg atoms are randomly placed in a single interac- 
tion volume as in Fig. 1(b), and the solution given by 
Eq. ([5]) is averaged over many random position realiza- 
tions (~ 10 4 ). 

The model interaction volume is chosen to be a cube of 
size -L, in order to assign random atom coordinates along 
the three Cartesian axes independently. This signifi- 
cantly simplifies the calculations, otherwise it is necessary 
to account for the possible correlations between the coor- 
dinates. Although actual excitation volume differs from 
a cube, special comparison between Monte Carlo simula- 
tions for the spherical and cubic volumes has shown that 
the results are almost identical if the excitation volumes 
are equal. Therefore, in the further calculations we use 
only cubic or parallelepiped volumes with homogeneously 
distributed atoms. 

For the correct comparison with Eq. (|10p we should 
properly define r$ for two Rydberg atoms in a cubic vol- 
ume V = I? . A natural way is to make this connection 
via the volume density of the two atoms in a cubic volume 
ri = 2/L 3 = (4 7 r^/3)- 1 that yields r « 0.5L. 

Figure 2 shows the numerical and analytical results 
obtained for the scalar (a) and full (b) DDI between two 
Rydberg atoms randomly placed in a cubic interaction 
volume. From Fig. 2(a) we see that usage of the scalar 
approximation results in the resonance amplitude that 
exceeds the saturated value 1/4 at #o ~ 2 — 14. The rea- 
son is a residual coherence due to incompletely damped 
Rabi oscillations, which are not taken into account in 
Eq. (fTDj) . Exact numerical solution of Eq. © reproduces 
well this peculiarity. This effect is completely absent in 
Fig. 2(b) where Monte Carlo simulation is done for the 
full DDI. Random orientation of the dipoles washes out 
the Rabi-like oscillations at any 6q, so that no coherence 
is expected for the Forster resonance between two disor- 
dered atoms. 
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FIG. 2: (Color online) Monte Carlo simulation of a Forster 
resonance between two Rydberg atoms randomly placed in a 
cubic interaction volume [as in Fig. 1(b)]. (a) and (b) De- 
pendence of the spatially averaged resonance amplitude on 
the average interaction pulse area 9q for the scalar (a) and 
full (b) dipole-dipole interaction; analytical fit has been done 
with Eq. (|1U|) . (c) and (d) Spectra of the Forster resonance 
in the dimensionless detuning scale A to at various 8o for the 
scalar (c) and full (d) dipole-dipole interaction. 



A very surprising finding is that Eq. (fT0|) fits well the 
numerical simulations in Fig. 2(b) at a = 0.55, although 
this equation has been derived for the scalar DDI. We 
suggest that Eq. (flOl) is a universal formula that can be 
used to obtain the spatially averaged Forster resonance 
amplitude in disordered atom ensembles. For example, 
a similar formula was used empirically in Ref. [14] to fit 
the time dependence of the number of interacting Ryd- 
berg atoms excited by a laser pulse. Universal scaling 
in a strongly interacting Rydberg gas was discussed in 
Ref. [27]. 

In Fig. 2 we also present the numerically calculated 
spectra of the Forster resonance in the dimensionless de- 
tuning scale A to for the scalar (c) and full (d) DDI ob- 
tained at various 0l that is, at various interaction vol- 
umes or interaction times. These spectra clearly show 
how the resonance saturates and broadens as 9q grows. 
Below saturation (f9 < 1) the ultimate full width at 
half maximum (FWHM) is equal to 27r, corresponding 
to 2n/to in the circular frequency scale. It is mainly a 
Fourier-transform limited width, which is larger than the 
average energy of DDI given by Og/to. For 8q > 1 the res- 
onance broadens and its amplitude saturates at the 0.25 
value. Averaging over the atom positions forms a reso- 
nance with broad wings and a cusp on the top, which is 
similar to that observed in atomic beams [28,29]. 

As we pointed out in Ref. [23], the resonance line 
shapes cannot be properly fit with a Lorentz profile. 
When we match the central part, the resonance wings 
are much broader than the Lorentz ones, due to rare in- 
teractions at very short distances. The commonly used 
Lorentz fit is thus inadequate for precise comparison be- 
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tween theory and experiment. 

Although the residuals of the Rabi-like population os- 
cillations remain in the spectra, from Fig. 2 we see that 
DDI in the disordered atom ensemble is incoherent and 
cannot be used to implement quantum phase gates. The 
two atoms must be spatially localized in the individual 
dipole traps to improve the coherence. 



D. Monte Carlo averaging for two localized atoms 

If the two atoms are localized in the individual dipole 
traps as in Figs. l(c)-l(d), the interaction strength de- 
pends on the spatial orientation of the traps with re- 
spect to each other and to the quantization axis z. The 
quantization axis is most convenient to choose along the 
dc electric field applied for Stark tuning of the Forster 
resonance and for ficld-ionization detection of Rydberg 
atoms. According to Eq. @, the strongest dipole-dipole 
interaction occurs when the two atoms are aligned along 
the z axis. In this case the axes of the two dipole traps 
in Fig. 1(c) are directed orthogonally to z, and spatial 
fluctuations occur mainly in the x — y plane. We are free 
to choose the dipole trap axes along the x direction for 
distinctness. In what follows we shall consider only this 
spatial configuration for the two atoms localized in the 
identical individual traps. 

Let us denote r x , r y , and r z to be the spatial uncer- 
tainties of the two atoms in the x, y, and z directions, 
correspondingly, and Rq to be the distance between the 
trap axes. The average interaction pulse area is defined 
for this configuration as 



_ 2V§tfcld23. 

7Q — ^— LQ. 

^ttsoHRq 
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We also need to choose the realistic spatial parameters 
of the dipole traps. Experiment [10] can be used for 
reference. The distance between the two traps was Rq rs 
10 /Lim. At the atom temperatures of ~ 200 /iK, the 
reported position probability distributions for each atom 
were approximately Gaussian with a ViZ ~ 0.3 /im and 
a x ~ 4 //m, corresponding to FWHM of r VlZ w 0.7 /zm 
and r x rs 9.4 /im. This dipole trap geometry can also be 
characterized by the aspect ratio r x /r VtZ w 13.4. 

Figures 3(a)-3(c) present our numerical Monte Carlo 
results obtained for the dipole trap geometry described 
previously, but the aspect ratio is varied from (a) 26.8 to 
(b) 13.4 and (c) 6.7. The right-hand pictures in Fig. 3 
show to scale the relative atom position uncertainties 
and the distance between the traps. Dependencies of 
the resonance amplitude on 8q (the left-hand panels in 
Fig. 3) demonstrate the damped Rabi-like population os- 
cillations. The first oscillation minimum around 9q « tt 
corresponds to the ir phase shift of the collective wave 
function. This point is what is needed to realize QPG. 
Nonzero minimum at this point results in a finite fidelity 
of QPG. Fidelity is given by the probability of the two 
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FIG. 3: (Color online) Monte Carlo simulation of a Forster 
resonance between two Rydberg atoms randomly placed in 
the two individual cigar-shaped [(a)-(c)] or symmetrical [(d)- 
(g)] dipole traps of various size. Left-hand panels: Depen- 
dence of the spatially averaged resonance amplitude on the 
average interaction pulse area 8o- Middle panels: Spectra of 
the Forster resonance in the dimensionless detuning scale A to 
for various 9q. The right-hand pictures show to scale the rela- 
tive atom position uncertainties and the distance between the 
traps. The z axis is oriented vertically. 



atoms to remain in the initial state 1 22) after this collec- 
tive state obtains the phase shift it: F = 1 — 2p2($o = 7r)- 
From Fig. 3 it is seen that fidelity (F 0.62) is 
poor for the trap geometry of Fig. 3(b) corresponding to 
the experimental conditions of Ref. [10]. An acceptable 
QPG fidelity of 0.85 can be obtained for the geometry 
of Fig. 3(c), that is, when the atoms are localized along 
the x axis two times better. This requires the atom tem- 
perature to be decreased from 200 /iK to ~ 50 fiK. The 
related spectra in Fig. 3 (the middle panels) show that 
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coherence is lost as Oq increases. One may conclude that 
the coherent dipole-dipole interaction of Rydberg atoms 
and its usage in QPG require the atoms to be colder and 
localized better than in the state-of-the-art dipole traps. 

We therefore have done the calculations for the sym- 
metrical dipole traps evenly localized in all directions 
[Figs. 3(d)-3(g)], which apparently can be realized at 
the lower atom temperatures. The highest coherence 
and fidelity of 0.96 are observed in Fig. 3(d) where 
r x,y,z = i?o/20. Such geometry corresponds to the 0.5 ^m 
in size dipole traps separated by 10 /jm. This is close 
to the ultimate spatial resolution achievable with opti- 
cal dipole traps. A somewhat worse fidelity of 0.86 is 
observed in Fig. 3(e) for the more realistic 1 /jm in size 
dipole traps. Larger traps of 2 /im (f) and 4 fxm. (g) in 
size demonstrate stronger decoherence (0.66 and 0.53 fi- 
delities, correspondingly) and are not suitable for QPG. 
Overall, Fig. 3 helps to estimate the attainable coherence 
and fidelity of QPG for various trap configurations. 



E. Monte Carlo averaging for dipole blockade 
between two atoms 



Dipole blockade (DB) must appear as the laser exci- 
tation of only one Rydberg atom out of a ground-state 
mesoscopic ensemble if the Rydberg atoms are strongly 
interacting [2]. It has been observed recently for just two 
Rydberg atoms interacting via quasi-Forster resonances 
in two close dipole traps [9,10]. 

In this section we will simulate DB at the Stark-tuned 
Forster resonance between two Rydberg atoms. There 
have been several articles on the numerical modeling of 
DB for many atoms (see, e.g., [30-37]). Most articles ana- 
lyzed the excitation probabilities in the resonance center, 
and only a few of them calculated the DB laser excitation 
spectra for some particular cases [31]. Here we present 
the extended calculations of the laser excitation spectra 
for two atoms in various spatial configurations, and an- 
alyze the effects of the Forster resonance detuning and 
laser non-monochromaticity. 

In order to calculate the dipole blockade spectrum for 
two atoms, we shall consider the narrowband laser ex- 
citation of Rydberg state |2) from ground state |0) [see 
Fig. 1(a)] for the two interacting atoms in the quasimolec- 
ular basis: 



* = ooo(*)|00) + oao(*)|20> + a a2 {t)\Q2)- 
022(t)|22)+oi3(t)|13)+a3i(t)|31). 



(12) 



Other collective states have large energy detunings and 
are not populated by the exciting laser radiation or DDI. 
The Hamiltonian is now taken as 



where E(t) is the electric field of the laser radiation. 
Laser-induced transitions occur from initial ground state 
1 00) to single-excited states 1 20) and |02), and to double- 
excited state 1 22). State 1 22) is coupled to other double- 
excited states 1 13) and |31), via DDI. Strong interaction 
shifts state 1 22) out of the resonance with laser radia- 
tion, and dipole blockade results in the decrease of the 
excitation probability of double-excited states. 

Probabilities to excite one or two atoms by a square 
laser pulse of duration i are calculated as 



Pi = ko(*o)r + K 2 (*o)n 

P2 = |a 22 (t )| 2 + |a 3 i(to)| 2 + |ai 3 (to)| S 



(14) 



H = 



(d a + 4) E{t) 



d a d b 
Ane 



1 



13l 

Ft 5 



(13) 



Blockade fidelity is given by P\. Perfect blockade corre- 
sponds to P\ = 1 and P2 = 0. 

In our Monte Carlo calculations the laser radiation is 
assumed to be linearly polarized along the z axis. Laser 
detuning 5 from the exact resonance |0) — > |2) is var- 
ied around its unperturbed zero value. Monte Carlo ex- 
citation spectra are calculated in the dimensionless de- 
tuning scale S tQ. The laser excitation pulse area is 
O = d()2Eoto/h, where ^02 is the dipole moment of transi- 
tion |0) -» 1 2) and Eq is the electric field amplitude of the 
laser radiation. If O = 7r and 5 = 0, the non-interacting 
atoms are completely transferred to Rydberg state |2). 
In the interacting atoms the probability of transfer to 
double-excited states reduces due to DB. 

In contrast to Eq. ([5]), there is no simple analytical 
solution to the Schrodinger equation for dipole blockade. 
Therefore, all further simulations have been done numeri- 
cally. Figure 4 shows the calculated Monte Carlo spectra 
of the one- and two-atom laser excitation probabilities 
Pi and P2 for O = ir and for various average interac- 
tion pulse areas 0q. Two Rydberg atoms are randomly 
placed in the two individual cigar-shaped [Figs. 4(a)-4(c)] 
or symmetrical [Figs. 4(d)-4(g)] dipole traps of various 
size or in a single excitation volume [Fig. 4(h)] [for the lat- 
ter configuration 8q is defined according to Eq. ©]. The 
right-hand pictures show to scale the relative atom po- 
sition uncertainties and the distance between the traps, 
which are the same as in Fig. 3. Forster resonance de- 
tuning A is taken to be zero. 

From the left-hand panels in Fig. 4 it is seen that there 
is almost no blockade at 9q = 2 for any trap configura- 
tion. These spectra are just probabilities to populate one- 
and two-atom states for the two non- interacting atoms, 
which are identical to what we observed experimentally 
at microwave transitions between sodium Rydberg states 
in our earlier work that relied on the atom-number re- 
solved study [15]. In the line center (5 = 0) probability 
to excite both atoms, P 2 , is close to unity, while proba- 
bility to excite one of the two atoms, Pi, has a dip and 
is nearly zero. Rabi oscillations appear in the resonance 
wings for our square laser excitation pulse. 

Dipole blockade becomes apparent at 80 = 8 (middle 
panels in Fig. 4). Probability P 2 reduces in the line cen- 
ter and there appears a dip corresponding to either par- 
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FIG. 4: (Color online) Monte Carlo simulation of the dipole 
blockade at a Forster resonance between two Rydberg atoms 
randomly placed in two individual cigar-shaped [(a)-(c)] or 
symmetrical [(d)-(g)] dipole traps of various sizes and in a 
single excitation volume (h). Spectra of the one-atom and 
two-atom laser excitation probabilities, Pi and P2, are shown 
in the dimensionless detuning scale 8 to for various average 
interaction pulse areas 9q. The right-hand pictures show to 
scale the relative atom position uncertainties and the distance 
between the traps. The z axis is oriented vertically. The 
laser excitation pulse area is O = it. The Forster resonance 
detuning A is zero. 



tial [(a), (b), (g), (h)] or almost complete [(c)- (f)] dipole 
blockade. At the same time, P\ grows and its dip con- 
verts to a peak. At 9q = 16 (right-hand panels in Fig. 4), 
these features transform to a strong blockade effect in 
the line center. In Figs. 4(c)-4(f) probability P2 drops to 
« 0.013 - 0.024, while Pi saturates at « 0.64. 

Large atom position uncertainties in Figs. 4(a) and 
4(h) substantially reduce the blockade efficiency at 9 n = 
16 due to zeros of the Forster interactions at certain 
dipole orientations [8] and due to lower interaction en- 
ergy between our z-dipoles in the x — y plane. Neverthe- 
less, we have found that blockade is still achievable for 



these spatial configurations. Blockade fidelity of 95% is 
reached at 9q > 100. This requires the interaction dis- 
tance to be R3 2 — 2.5 times smaller than that for the well 
localized atoms in Figs. 4(c)-4(f). 

Further analysis is done with Fig. 5. The left-hand 
panels in Fig. 5 demonstrate that the Pi value can be 
increased up to 1 if we set a smaller laser pulse area 
G = This reduction is due to the yf~N depen- 

dence of the collective Rabi frequency on the number of 
interacting atoms N [30]. In Fig. 5 we use this new 
value in all simulations. 

The middle panels in Fig. 5 show the dependencies 
of the blockade fidelity (Pi value) on the average dipole 
interaction pulse area 80 for three different detunings 
from the exact Forster resonance: A = 0, 207r/to, and 
407r/i . These correspond to A/2tt =0, 10, and 20 MHz 
detunings at the to — 1 fJS interaction time. Zero de- 
tuning is relevant to the precisely Stark-tuned Forster 
resonances, while 5-20 MHz detunings are typical for 
accidental quasi-Forster resonances [8]. It is seen that 
dipole blockade has a threshold-like behavior versus 9o, 
and non-zero detuning substantially increases the thresh- 
old. Nevertheless, the detuning effect can be overcome 
by a stronger interaction, so that it is still possible to 
reach Pi « 1. This has been confirmed experimentally 
in Refs. [9,10]. 

Figures 4 and 5 help us to understand what kind of 
one- and two-atom excitation spectra can be observed in 
the presence of DB. We note, however, that these calcula- 
tions were made for a monochromatic laser radiation (the 
laser linewidth Aujl *C 2n/to) and square laser excita- 
tion pulse, which produce the Rabi-like oscillations in the 
spectra. The oscillations are washed out if Aujl increases 
due to technical imperfections or if parasitic electromag- 
netic fields broaden the resonance. This leads to decoher- 
ence and decreases the blockade fidelity. The right-hand 
panels in Fig. 5 show the dependencies of Pi on 6q for 
the three different laser lincwidths: Aljl = 0, 2ir/to, and 
47r/£o- These correspond to Aojl/^tt =0, 1, and 2 MHz 
at the to = 1 /is interaction time. Finite laser linewidth 
in our Monte Carlo simulation was introduced as ran- 
dom laser frequency jumps around the zero value. It is 
seen that laser line width inevitably decreases the dipole 
blockade fidelity, which cannot be improved by a stronger 
interaction. Therefore, depending on the linewidth of the 
available laser, it is crucial to use to <C l-nj Alol in the 
experiments in order to achieve high blockade and quan- 
tum gate fidelities. 



III. INTERACTION OF MORE THAN TWO 
RYDBERG ATOMS 

A. Frozen atoms 

Let us now consider a Forster resonance between N > 
2 frozen Rydberg atoms with |1), |2), and |3) Rydberg 
levels participating in the resonance [see Fig. 1(a)]. The 
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FIG. 5: (Color online) Monte Carlo simulation of the dipole 
blockade at a Forster resonance between two Rydberg atoms 
randomly placed in two individual cigar-shaped [(a)-(c)] or 
symmetrical [(d)-(g)] dipole traps of various sizes and in a sin- 
gle excitation volume (h). The left-hand panels demonstrate 
the dipole blockade spectra at the laser pulse area 
The middle panels show the dependencies of the blockade fi- 
delity Pi on the average dipole interaction pulse area 80 for 
three different detunings from the exact Forster resonance: 
A = 0, 20-K/to, and 407r/to- The right-hand panels show the 
dependencies of Pi on 80 for three different laser linewidths: 
Awl = 0, 2n/to, and 4tt /to- The right-hand pictures show to 
scale the relative atom position uncertainties and the distance 
between the traps. 

main calculation difficulty is that basis size Z(N) of the 
collective wave function rapidly grows with N. It is equal 
to the Gegenbauer polynomial C^ N \l/2), which yields 
Z{2) = 3, Z(3) = 7, Z(4) = 19, Z(5) = 51, Z(6) = 141, 
... , Z(10) = 8953, ... . 

We assume that at t = a short laser pulse ini- 
tially excites the collective state |22222...2). This state 
is coupled via Forster resonance to Nl/(N — 2)1 states 
of the kind 1 13222. . .2). In turn, these states are cou- 
pled to Nl/[(N - 4)!(2!) 2 ] states of the kind |13132...2), 



and so on. Taking into account all possible permuta- 
tions, a general formula for the number of such states 
is n k = N\/[(N ~ fc)!(0/2)!) 2 ], where k is the number 
of Rydberg atoms that have abandoned state |2) due to 
DDI (k is always even for the Forster resonance, because 
one atom goes to the lower Rydberg state and the other 
atom simultaneously goes to the upper state) . The tran- 
sition probability per atom, or fraction of atoms in final 
state |1), is calculated as 

k ™ fc 

p n= ^y5Ii a * fc i 2 ' ( 15 ) 

fc=2,4,6... i=l 

where a$ & are the probability amplitudes of collective 
states, and k/(2N) is a normalization factor for each k. 
The maximum k value is equal to N for the even N, or 
to N - 1 for the odd N. 

It is impossible to derive a simple analytical for- 
mula like Eq. ([5]) for pn at N > 2, even for frozen 
Rydberg atoms. This is due to the fact that be- 
sides the resonant dipole-dipole interaction Rb(37P) + 
Rb(37P) -> Rb(37S) + Rb(385), it is necessary to take 
into account the exchange processes like Rb(37P) + 
Rb(37S) -> Rb(37S) + Rb(37P) and Rb(37P) + 
Rb(38S) Rb(385") + Rb(37P), which are always reso- 
nant (independent of the electric field) and lead to a pop- 
ulation hopping between neighboring atoms [17]. This 
hopping broadens and decoheres the Forster resonances, 
although it is not the main source of the broadening [38] . 

If, however, we consider a weak dipole-dipole interac- 
tion, most of the population remains in the initial state 
|22222...2), and only states with k = 2 are weakly popu- 
lated. In this case the always-resonant processes can be 
neglected, since S states in each atom are weakly popu- 
lated too. Then it can be shown that Eq. ([5]) is simply 
modified by replacing 2£l 2 ab with the sum of the squares 
of the interaction energies for all atomic pairs: 

and Eq. (TT5|) finally yields 

^^r^^l^^'")- (17) 

Exact numerical modeling with the exchange interac- 
tions fully taken into account has confirmed that Eq. (JTTJ) 
works well at pat (to) < 0.1, both in the resonance center 
and wings. 

For the strong DDI the maximum pn value depends on 
the N evenness. It is 1 /2 for the even N and (JV— 1)/(2JV) 
for the odd N. In particular, for three atoms the maxi- 
mum transfer probability is 1/3, since only one of the 
three atoms can be transferred to state |1) at a time. De- 
coherence due to spatial fluctuations may decrease these 
maximum attainable values by up to 2 times. 
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B. Spatially fluctuating atoms 

We can imagine various spatial configurations includ- 
ing linear, planar, cubic, honeycomb, or ring lattices of 
single ground-state atoms intended for a quantum com- 
putation. Corresponding considerations were carried out 
in many theoretical papers (see, e.g., [30,32,33]). 

Two-qubit logic gates imply simultaneous laser excita- 
tion of only two Rydberg atoms. Our results obtained in 
the previous section can be applied for their analysis with 
minor changes concerning mainly the renormalization of 
the interaction strength. 

In this section we consider simultaneous laser excita- 
tion and interaction of more than two Rydberg atoms. 
One may expect that many-body phenomena would en- 
hance the DDI strength in the mesoscopic ensembles, but 
coherence may be lost upon averaging over many inter- 
acting atoms. 

The numerical Monte Carlo simulations were per- 
formed for iV = 2 — 5 interacting Rydberg atoms, initially 
excited to state |2) and randomly positioned in a single 
cubic or cigar-shaped volume or in the evenly spaced indi- 
vidual dipole traps (see the right-hand pictures in Fig. 6). 
In this approach, the time evolution of all Z(N) collective 
states is obtained by numerically solving the Schrodinger 
equation for various N, A, and interaction times t$. We 
accounted for all possible binary resonant interactions be- 
tween N atoms, as well as the always-resonant exchange 
interactions. The initial positions of N atoms were aver- 
aged over 500 random realizations. A similar approach 
was used in Refs. [21,22]. 

Figure 6(a) shows the results of our calculations of the 
Forster resonance amplitude and spectrum for two to five 
Rydberg atoms randomly positioned in a cubic volume. 
The dependencies of the amplitude on the average inter- 
action pulse area 0q [defined for two atoms according to 
Eq. ©] demonstrate similar saturation curves for all N. 
The curves for N=2, 4, and 5 saturate at pn — 0.25 and 
are fitted well with Eq. (TTU1) at a 2 ~ 0.55, ps 1.60, and 
a 5 w 2.10. 

The case N = 3 is distinguished by a different satura- 
tion amplitude P3 ss 0.21 and a different saturation curve, 
which is similar to that in Fig. 2(a). The saturated am- 
plitude lies between the maximum possible value of 1/3 
for coherent interaction and 1/6 for completely incoher- 
ent interaction. This indicates that residual coherence 
is always present in the interaction of three atoms. In 
general, this observation should be typical for the odd 
number of the interacting Rydberg atoms. For example, 
the resonance amplitude for five atoms in Fig. 6(a) sat- 
urates at p5 fa 0.25 and also lies between the maximum 
possible value of 2/5 for coherent interaction and 1/5 for 
incoherent interaction. 

The spectra in Fig. 6(a) are calculated for the unsat- 
urated (0o = 0.5) and saturated (9q = 8) Forster reso- 
nances. In the unsaturated spectra the ultimate FWHM 
is defined by the inverse interaction time 2ir/to for all 
N. It is mainly a Fourier-transform-limited width, unless 




FIG. 6: (Color online) Monte Carlo simulation of the Forster 
resonance between N = 2 — 5 Rydberg atoms randomly placed 
in (a) a single cubic excitation volume, (b) a single cigar- 
shaped dipole trap, and (c) individual evenly spaced sym- 
metrical dipole traps. Left-hand panels: Dependence of the 
spatially averaged resonance amplitude on the average inter- 
action pulse area between two atoms 9q. Middle panels: Spec- 
tra of the Forster resonance in the dimensionless detuning 
scale A to at 9o = 0.5 and 8o = 8. The right-hand pictures 
show to scale the relative atom position uncertainties. The z 
axis is oriented vertically. 



&N0o < 1< The probability for two atoms to interact 
at shorter distances (where a^0Q > 1) is low, and the 
resonance is narrow despite the spatial averaging. The 
signatures of the Rabi-like oscillations are visible in the 
resonance wings. The resonances broaden as N grows, 
due to the increase in the average energy of DDI. 

Figure 6(b) presents the amplitude and spectra cal- 
culated for two to five atoms confined in a single cigar- 
shaped dipole trap oriented along the z axis. The as- 
pect ratio of the trap is r z /r XtV = 10 (for example, a 
10-/im-long dipole trap with l-//m transverse size, as in 
Ref. [10]). As it was pointed out in Ref. [39], mainly 
nearest-neighbor atoms interact in this specific linear ge- 
ometry. This reduces the saturated resonance amplitude 
to a value below 0.25, because for TV atoms in the elon- 
gated dipole trap there are less than N effectively inter- 
acting atoms on the average. 

Figure 6(c) shows the amplitude and spectra calcu- 
lated for the linear dipole trap arrays of N = 2 — 5 atoms 
aligned along the z axis. The dipole traps are evenly 
spaced at distance r, while the atom position uncertain- 
ties are r/5 [for two atoms this geometry corresponds to 
Fig. 3(f)]. The value of 0q is defined for two neighboring 
atoms according to Eq. (jlip. As the atoms are local- 
ized well in all directions, the interaction is coherent to 
some extent. The maximum coherence is observed for 
two atoms. At 0q = 0.5 the Forster resonance spectra 
are almost identical for all N due to primarily nearest- 
neighbor interaction. At 0q = 8 there appear the asym- 
metries in the spectra which have different characters for 
different N. 
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FIG. 7: (Color online) Monte Carlo simulation of the dipole 
blockade at a Forster resonance between three Rydberg atoms 
randomly placed in a single cubic excitation volume. Spectra 
of the one-atom, two-atom and three-atom laser excitation 
probabilities Pi, P2, and P3 are shown in the dimensionless 
detuning scale S to for various average interaction pulse areas 
6*o. Laser excitation pulse area is G = tt [(a)-(c)] and & = 
7r/\/3 (d). The Forster resonance detuning A is zero. 



Figure 6 demonstrates that in all spatial configura- 
tions coherence is lost as N increases, although the DDI 
strength substantially grows. We therefore cannot bene- 
fit from increasing the number of the interacting Rydberg 
atoms in performing the quantum phase gates. 

On the other hand, dipole blockade is less sensitive 
to decohcrcncc and it might benefit from increasing N. 
With our computer code we are able to calculate the DB 
excitation spectra and fidelity for up to five atoms in 
various spatial configurations. 

As an example, Fig. 7 presents the Monte Carlo results 
for the dipole blockade at a Forster resonance between 
three Rydberg atoms randomly placed in a single cubic 
excitation volume. Spectra of the one-, two-, and three- 
atom laser excitation probabilities P\, P2, and P3 are 
shown in the dimensionless detuning scale S to for various 
average interaction pulse areas 9q [defined for two atoms 
according to Eq. ©]. The laser excitation pulse area is 
6 = tt in Figs. 7(a)-7(c) or 6 = in Fig. 7(d). The 

Forster resonance detuning A is zero. 

From Fig. 7 we see that the probability to excite three 
atoms, P3, rapidly drops as 9q increases, but probability 
to excite two atoms, P 2l drops much slower. This means 
that for the same atom density the DB does not bene- 
fit much from increasing N. Blockade becomes stronger 
when we place more atoms into the same interaction vol- 
ume, that is, when we increase the atom density in the 
mesoscopic ensemble. Another theoretical study [32] has 
also revealed that at N = 3 dipole blockade cannot be 
complete for certain spatial configurations. 



ade between two to five cold Rydberg atoms. Various 
spatial configurations have been analyzed in the context 
of quantum-information processing with neutral atoms in 
optical dipole traps. 

It has been shown that quantum phase gates can be 
implemented using coherent Rabi-like population oscilla- 
tions at exact Forster resonance. However, this can be 
done only with well localized atoms at the temperatures 
below 50 [iK. Higher temperatures lead to larger atom 
position uncertainties and to decoherence of the popula- 
tion oscillations. Increasing the number of the interacting 
Rydberg atoms also leads to decoherence and does not 
help in implementing the quantum phase gates. 

Dipole blockade at the Stark-tuned Forster resonance 
between few Rydberg atoms is more robust against the 
position uncertainties. The dependence of the blockade 
fidelity on the dipole-dipole interaction strength demon- 
strates a threshold-like behavior. The threshold value 
substantially increases if the Forster resonance detuning 
is not zero, but blockade is still achievable at the larger 
interaction strengths. In contrast, finite line width of 
the exciting laser radiation decreases the blockade effi- 
ciency, which cannot be improved by a stronger interac- 
tion. Increasing the number of Rydberg atoms does not 
increase blockade fidelity, since for the same atom density 
the two-atom excitation probability is almost the same 
as for just two atoms in the interaction volume. Finally, 
we have shown that dipole blockade can be observed even 
with two disordered atoms in a single excitation volume, 
if the interaction strength is large enough. 

We conclude that both quantum phase gates and 
dipole blockade have their advantages and drawbacks 
and they require different experimental conditions to be 
observed. In particular, quantum phase gates can be 
performed at a weak dipole-dipole interaction, but the 
atoms should be localized well to preserve coherence. 
Dipole blockade is insensitive to position uncertainties, 
but suffers from the finite laser linewidth and requires 
much stronger interaction, especially for quasi-Forster 
resonances. The results obtained can be useful for the 
analysis of the experimental conditions appropriate to 
quantum-information processing with Rydberg atoms. 
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